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CTN ' The unusual slow carrier relaxation dynamics - observed in femtosecond pump-probe experiments 

on high-temperature superconductors and recently also in a charge-density-wave system - is analyzed 
in terms of a model for relaxation of carriers in intra-gap states. The data on YBa2Cu307_i near 
optimum doping and K0.3M0O3 are found to be described very well with the model using a BCS-like 
gap which closes at Tc- From the analysis of the data we conclude that a significant intra-gap density 
of localized states exists in these materials, which can be clearly distinguished from quasiparticle 
states by the time-resolved optical experiments because of the different time- and temperature- 
dependences of the photoinduced transmission or reflection. Localized charges are suggested to be 
the most likely origin of the intra-gap states, while the similarity of the response in the two materials 
appear to exclude spin and vortex excitations. 
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I. INTRODUCTION 

Photoinduced absorption or reflection spectroscopy us- 
ing femtosecond lasers is potentiaUy a very powerful tool 
for the study of the electronic structure of supercon- 
ductors and related materials. Recent pump-probe ex- 
periments on cuprate superconductors [|]-Q and some 
other materials have shown that a photoinduced 
change in absorption or reflection can be observed at 
low temperatures and especially for T < Tc- The ef- 
fects are believed to be caused by excited state absorp- 
tion of the probe pulse from photoexcited quasiparticlc 
(QP) states B and theoretical analysis of the response 
was found to be in good agreement with experimen- 
tal data on YBa2Cu3 07_5 over a wide range of dop- 
ing W. However, in addition to the QP response which 
occurs on the picosecond and subpicosecond timescale, 
a distinct slower response was also consistently ob- 
served in YBa2Cu307-5 (YBCO) J|, BiaSrzCaCuaOg 
(BISCO) and Bi2Y2;Cai_a;SrCu208 Q and more recently 
in the charge-density- wave (CDW) quasi one-dimensional 
Peierls insulator K0.3M0O3 0. It was thought to be of 
non-thermal origin (detailed analysis is given in Ref. g]) 
and occurs on a timescale of 10~^ s or longer (see Fig. 
1 for details). Its anomalous T-dependence [|| - which 
is qualitatively different from the T-dependence of the 
fast QP recombination dynamics - lead the authors to 
the suggestion that the signal is due to localized states 
near the Fermi energy. However, the processes involved 
were so far not discussed in any detail. 

In this paper we examine quantitatively the photoin- 
duced absorption (or reflection) from localized intra-gap 
states. We develop a theoretical model for two different 
cases: (i) the case of a material with a BCS-like collec- 
tive gap and (ii) the case of a T-independent gap (some- 
times called a pseudogap) which exists above Tc- We 
compare the model predictions with the available data 
on the temperature-dependence of the photoinduced re- 
laxation on the nanosecond timescale in two different 
materials, both of which are generally thought to have 
a low-temperature gap: YBa2Cu306.9 {Tc = 90K) and 
K0.3M0O3 (T, = 183K). 



II. THEORETICAL MODEL 

As the name implies, pump-probe spectroscopy in- 
volves the excitation of the material by an ultrashort 
pump laser pulse and the subsequent measurement of 
the resulting change in optical absorption, transmission 
or reflection of the sample caused by photoexcited charge 
carriers. As photons from the pump laser pulse are ab- 
sorbed, they excite electrons and holes in the material 
(see schematic diagram in Fig. 2 a)). These particles re- 
lease their extra kinetic energy by scattering amongst 



themselves and with phonons (step 2 in Fig. 2a)). This 
energy relaxation process is very rapid and the particles 
end up in QP states near the Fermi energy {Ep) within 
Te =10~100 fs Q. Subsequent relaxation is slowed down 
by the presence of the gap and a relaxation bottleneck 
is formed. From pump-probe photoinduced transmission 
experiments in YBCO |i|,| and BISCO §, the relax- 
ation times of the photoexcited QPs were found to be in 
the range tqp =0.3~3 ps. 

In addition to the picosecond transient, the signal on 
nanosecond timescale has been consistently observed in 
HTSC [|J6|jl| and recently also in quasi ID CDW insula- 
tor K0.3M0O3 [g. In Fig. 1 the photoinduced transient 
taken on YBCO at T=80 K is shown (squares). After 
photoexcitation (i = ) the signal relaxes within 10 ps 
to some non-zero value, that can be represented by the 
constant on the 100 ps timescale |6). The lifetime of the 
slow component, r^, cannot be directly measured, since 
it appears to be longer than the inter-pulse separation of 
tr ~ 10 — 12 ns. This results in the signal pile-up due 
to accumulation of the response from many pulses. The 
magnitude of this pile-up, j4, given by the difference in 
the signal amplitude when the pump pulse is unblocked 
(2), and zero signal, when the pump is blocked (1) at 
negative time delays (see Fig. 1), is several times larger 
than single pulse contribution, ag, given be the difference 
between (4) and (2). 

Assuming - for simplicity - that relaxation of the signal 
is exponential, we can write an equation for the steady 
state amplitude A: 

00 
A = ao ^ exp{-ntr/TL) = ao/(l - exp{-tr/TL)) (1) 

n=0 

Here oq is the amplitude of the signal from a single pulse 
and tl is the relaxation time. Usually, A ^ ao and we 
can expand the exponent in the denominator of Ecl(1), 
to obtain tl — tr{A/ao). From the experiments |g,^ it 
appears that tl > 10^^ s. This is long in comparison 
with the phonon relaxation time and with the phonon 
escape time from the excitation volume into the bulk or 
thin film substrate, which is typically lO^^'^s, so we can 
ignore phonon escape effects and discuss only intrinsic 
relaxation processes. 

The process giving rise to the actual photoinduced op- 
tical signal from intra-gap states is shown as step 3 in 
Fig. 2a). Here we do not discuss the optical probe pro- 
cess in detail, but make the very general assumption that 
the photoinduced change in sample transmission A T / T 
or reflection ATZ/TZ is proportional to the photoinduced 
density of localized states populated by the laser pump 
pulse. The photoinduced signal AT/T (or ATZ/TZ) is 
then proportional to the number of filled localized states 
and the T-dependence is mainly determined by the oc- 
cupation of the intra-gap states. 

Since the relaxation time tl is long in comparison 



with tqp and phonon relaxation times, we assume that 
phonons and quasiparticles can be described by equihb- 
rium densities N^ and N respectively. For the relaxation 
of the localized carriers we apply arguments similar to 
those originally proposed by Rothwarf and Taylor 0. 
The rate equation for the total density of localized exci- 
tations Nl is then given by: 



dt 



-RNl - ^Nl +-fN + f3N^ 



(2) 



The first term in Eq.(2) describes the recombination of 
two localized excitations to a Cooper pair with a recom- 
bination rate R. The second and the third terms de- 
scribe the exchange of an electron or a hole between the 
localized and quasiparticle states with density N^ and 
N respectively and with a rates 7/7 ex exp {— AE / k gT) 
where AE is the energy barrier between trapped carriers 
and quasiparticles |^. The last term describes the spon- 
taneous creation of localized excitations by phonons with 
a relaxation rate /3. The four processes in Eq. (2) are 
shown schematically in Fig. 2b)). 

Assuming the ansatz for iV^ = A^^o + "^l , where iV^g 
is equilibrium density of localized particles and nj^ is the 
photoinduced density created by the laser pulse and tak- 
ing into account that iV^ and N are given by their equi- 
librium values, we can rewrite Eq.(2): 



^ - -Rnl - {2RNlo + l)nL 
at 



(3) 



This equation is sufficiently general that it can be applied 
to different superconductors and different gaps. (For the 
case of a CDW gap, the Cooper pairs are replaced hy e — h 
pairs.) In this paper we consider a) a superconductor 
with a BCS-like gap (which can also be used in the case 
of a gap formed by a CDW) and b) a Bose condensate 
of pre-formed pairs with a T-independent gap. In the 
latter case, the gap is better considered as an energy- 
level splitting between paired and unpaired states, and 
the main difference is that the gap is T-independent and 
exists above Tc. A general exact solution to Eq.(3) is 
given in the Appendix. Next we consider the limiting 
behavior at low T and near Tc- 

In the case of a collective BCS-like gap, A^ (T) , the 
recombination rate below Tc is to the lowest order in A^ 
proportional to the square of the order parameter P,10| : 



R 



,{A,{T)/ncf 



(4) 



where a is a constant, and f2c is the phonon spectrum 
cutoff frequency. (Above Tc all recombination processes 
disappear as As{T) -^ 0.) To relax the carriers in local- 
ized intra-gap states via quasiparticle states (the term 
jriL in Eq.(3)), an energy of the order of A^ is re- 
quired and so this process is exponentially suppressed 
at low temperatures. The term proportional to A^^o 



is also small at low temperatures, because the num- 
ber of thermally excited localized excitations is small as 
{kBT/As{0))'^, where fi depends on the density of local- 
ized states. Therefore, for analysis of the relaxation of 
localized excitations at low temperatures T <C Tc, we re- 
tain only the first term in Eq.(3) giving a solution of the 
form: 



nLit)=nL{0)/{nLiO)Rt + l) 



(5) 



To obtain the stationary solution for a repetitive laser 
pump pulse train excitation, we use the condition that 
the total number of localized excitations that recombine 
between two laser pulses should equal the number of lo- 
calized excitations created by each laser pulse: 



riLiO) ~ nLiO)/inLiO)Rtr + 1) = mphiT) 



(6) 



where 77 oc ^tqp — if /A (T) is the probability of trapping 
a QP into a localized state and nQp(T) is the number 
of photoinduced QPs at temperature T created by each 
laser pulse. Since the number of photoexcited carriers 
is typically small compared to the overall carrier density 
'nnqp <^ ni(0), we can estimate ^^(0) as 



riLiO) 



lr]nQp{T) 



Rtr 



(7) 



As a result, combining (4) and (7) we get an expression 
for the T-dependence of the photoinduced transmission 
(or reflection) amplitude at low temperatures: 



|Z\T/r|ocnL(0) 



lr]nQp{T) ric 



atr A,(T) 
where ngp for a BCS-like case is given by M: 

.rr.. £,/(A,(T)+fcBT/2) 

nqpiT) = — 
1 



N{0)M1a 



-J^^eM-^s{T)/kBT) 



(8) 



(9) 



Here Sj is the energy density deposited per pulse, v is 
the effective number of phonons per unit cell involved in 
the relaxation process and N{0) is the density of states 
at the Fermi energy in units eV~^cell~^spin^^. 

Near Tc, when the number of thermally excited local- 
ized carriers becomes comparable or larger than num- 
ber of nonequilibrium carriers, the relaxation terms 
2RNLonL and "/ul become dominant, and the solution 
to Eq. (3) has an exponential rather than a power law 
time dependence of the form nL{t) — 71^(0) exp {—t/r) , 
where 1/r = {2R{T)Nlo + 7) . For T ^ T^, the temper- 
ature dependence of the photoinduced signal amplitude 
is then given by: 



|ZiT/T|(X7ii(0) = 



WQp{T) 



{2NloR{T) + 7)i, 



(10) 



Thus the predicted amphtude of the signal increases with 
increasing T up to Tc and drops to zero above Tc- 

In the case of a T-independent gap (or "pseudogap"), 
AP, we assume that the gap exists at all temperatures. 
Since the gap does not close at Tc, the recombination 
rate does not go to zero at Tc and instead of (4) we have 
a constant, F. 

However, below Tc the presence of the condensate may 
also have an effect on the recombination of localized ex- 
citations. In general, the relaxation rate is a function of 
the order parameter. To take this into account, we can 
expand it in terms of even powers of A. Near Tc the or- 
der parameter is small and we can keep only the lowest 
power in A^ [[l3[ . If we assume that the order parameter 
exhibits mean-field behavior (A oc ^Jl — T/Tc), then: 



i?~a(l-T/Tc)-t-r 



(11) 



where a is phenomenological constant which describes 
the dependence of the relaxation rate on T below Tc, 
which in general is not equal to 0. 

A slightly different expression for R is obtained if we 
assume that the recombination rate is dependent on the 
pair momentum in the condensate through the kinetic 
energy. To illustrate this, we write the recombination 
rate as F = Fq -|- F where F is the momentum averaged 
recombination rate, and Fq is the recombination rate for 
pairs with fc = 0. The recombination rate is generally 
proportional to the number of pairs in the condensate 
Up. For Bose condensation this is given by rip oc (1 — 
{T /TcY^'^), and we thus obtain a formula for the total 
relaxation rate which is similar to Eq. (11), but with a 
different temperature dependence in the first term: 



i? = (Fo-f)(i-(r/rc)3/2) + f 



(12) 



In principle the two cases (11) or (12) can be distin- 
guished by measurements of the T-dependence of the 
photoinduced transmission or reflection below Tc, al- 
though the difference will be very small and very high 
quality data is needed to do this. 

To obtain the photoinduced signal amplitude, we sub- 
stitute Eq.(ll) into eq. (3): 



\AT/T\ (xni(O) 



''QP 



{a{l-{T/TcY) + T)U 



-, (13) 



where /3 is a constant (generally 1 or 3/2) and Upp is the 
equivalent of formula (9) for the case of a temperature 
independent gap W: 



^QP 



Ei/AP 



1 



2i/ 



N{o)nn, 



■ exp(- AP/fcsT) 



(14) 



In contrast to the BCS case, expression (13) is non-zero 
above Tr and reduces to: 



|zir/r|(xni(o) = 



"-QP 



Ttr 



(15) 



which implies that the photoinduced absorption signal 
should remain observable well above Tc and should reveal 
the presence of a pseudogap if it exists. 

Just as before, a crossover to exponential time relax- 
ation takes place when the number of thermally excited 
excitations become large and the second term in the 
Eq.(3) becomes dominant, i.e. T ~ A^, leading to a lin- 
ear intensity dependence of \AT/T\ given by modified 
Eq.(lO) with Uqp and R from Eq.(12). Note that this 
crossover may occur at higher temperatures than exper- 
imentally measured since T-independent gap is typically 
AP > 300K. 



III. COMPARISON WITH EXPERIMENTAL 
DATA 

The temperature dependence of \AT/T\ or \ATZ/TZ\ 
for optimally doped YBa^CusOy-a [ijlll and for 
K0.3M0O3 g] is plotted in Fig. 3a) and b) respectively. 
In both compounds the signal amplitude increases with 
increasing temperature followed by an abrupt drop above 
Tc- Now we compare the predicted T-dependences of 
\AT/T\ (Eqs. (8) and (10)) with the data. The param- 
eters used in the fits are the same as previously used 
in the analysis of the fast relaxation component [Q|^] 
with values of dimensionless constant jy/^-^^ ~ 30 for 
YBaaCugOr-a @ and ~ 10 for K0.3M0O3 ||!' Since the 
magnitude of the gap is of the order of Ac (0) ~ 5fcTc we 
expect that Eq. (8) is valid up to temperatures close to 
Tc- In Fig. 3a) and b) we show the calculated tempera- 
ture dependences of \AT/T\ using Eq.(8) in comparison 
with experimental data for YBa2Cu307_5 and K0.3M0O3 
respectively. We would like to stress that Eq.(8) is inde- 
pendent of the shape of DOS of localized states within 
the gap and shows a universal temperature dependence. 
It can be seen from the Fig. 3b) that at low temperature 
there is deviation of the calculated curve from the experi- 
mental points. To explain this effect we should remember 
that Eq. (4) for constant R is valid near Tc where Ac (T) 
is small. (If we add next - fourth order - term in the 
expansion of R in powers of Ac (T) we can account for 
this discrepancy.) 

Near Tc Eq.(8) fails and Eq. (10) should be used (see 
also Appendix) . It leads to a crossover from square root 
intensity dependence of \AT/T\ at low T described by 
Eq.(8) to linear intensity dependence near Tc predicted 
by Eq.(lO). 

Finally let us discuss the effect of 7. In Fig. 3a) and b) 
fits to the data using the general solution (Eq.(19)) are 
shown with dashed lines. In these fits we have also added 
the fourth order term in the expansion of R in powers of 
Ac (T). As can be seen from these fits the effect of 7 



becomes important near Tc by cutting the divergence of 
\AT/T\ as T^Tc. 

In Fig. 3c) we have plotted calculated values of 
\AT/T\ oc nL{0) for the case of temperature indepen- 
dent pseudogap (which might be applicable in under- 
doped cuprates for example) as a function of temperature 
for different values of parameter a using Eq.(13). As can 
be seen from the Fig. 3c) in this case slow relaxation via 
localized states is present also above Tc- This effect is 
due to temperature dependence of Jigp controlled by T- 
independent pseudogap above Tc. 

In cuprate superconductors there is spectroscopic ev- 
idence suggesting that there is a significant density of 
states in the gap possibly extending to the Fermi level, 
which is often attributed to a d-wave gap symmetry. 
However, by normal spectroscopies it is difficult to de- 
termine if the states in the gap are QP states or, for 
example, localized states. Time-resolved spectroscopy 
can answer this question rather effectively because of the 
different time- and temperature-dependences of the QPs 
and localized carrier relaxations. It was argued that in 
the presence of impurity scattering QP DOS in the d- 
wave state remains finite at zero energy |I3. Recently 
it was proposed [ [L5[ that the quasiparticles in the super- 
conducting state may become strongly localized for short 
coherence length d-wave superconductors. However this 
statement has been questioned by Balatsky and Salkola 
pM and remains controversial. On the basis of available 
experimental data we cannot make any definite conclu- 
sion about origin of intragap localized states. 

We can, however estimate the density of the intra-gap 
states from the available data by assuming that the opti- 
cal probe process (step 3) is similar for excited state ab- 
sorption from localized states and for QPs. Both optical 
probe processes involve allowed transitions to the same 
final state E2 and so this assumption is not unreasonable. 
From typical photoinduced reflection data for YBCO (as 
in Figure 1), we find that approximately \AT/T\j^ ~ 
\AT/ T\qp, implying that also 71l(0) ~ nqp. From this 
we can conclude that the density of intra-gap states is 
comparable with the density of QP states. This obser- 
vation has important implications for the interpretation 
of frequency-domain spectroscopies, since it suggests that 
the spectra should show a very significant intra-gap spec- 
tral density due to localized states, irrespective of the gap 
symmetry. 

Assuming that the optical transition probability of the 
probe pulse is the same for QPs as for the intra-gap 
states, from Eqs. (8) and (9) we obtain: 



\AT/T\ 

\atJt\. 



riL 
nqp 






(16) 



where \AT / T\ is the photoinduced transmission due to 
the QPs. From Fig. (1) typically riL/nqp = 0.1 — 1 and 
using a pulse repetition rate tr = 12 ns and with tl — 100 



ns, we obtain an estimate of the trapping probability for 
carriers by localized states of 77 = 0.1 — 1. 

A detailed discussion of the origin of the localized 
intra-gap states in the cuprates should be deferred until 
more systematic data as a function of doping is available, 
and we only mention some of the most likely possibilities: 
(i) localized states associated with the inhomogeneous 
ground state of the cuprates (stripes) ||l^ , (n) intrinsic 
defect states, (iii) localized QP states in d-wave supercon- 
ductor [^ and possibly (iv) holons |18|. In K0.3M0O3, 
the nature intra-gap excitations has been a subject of ex- 
tensive study over the years and the reader is referred to 
ref. 1 12 1 for a review. However, the fact that the signals 
in K0.3M0O3 and YBa2Cu307_5 are very similar appears 
to rule out both spin excitations and vortex states, leav- 
ing localized charges as the most plausible origin of the 
intra-gap states. 



IV. CONCLUSIONS 

To conclude, the calculated time- and temperature- de- 
pendence of the photoinduced absorption for the case of 
a BCS-like gap is found to be in good agreement with 
experimental data from femtosecond time-resolved spec- 
troscopy on the cuprate superconductor YBa2Cu307_5 
near optimum doping and the charge-density wave insu- 
lator K0.3M0O3. We find that time-resolved spectroscopy 
can very effectively distinguish between QP states and 
localized states in the gap. A rather surprising feature 
of the data is the remarkable separation of the QP re- 
sponse on the femtosecond timescale and the slow re- 
sponse of intra-gap state relaxation on the scale of 100s 
of nanoseconds. In both materials we find a significant 
intra-gap density of states, which display very different 
time-dynanfics and T-dependence than the QP states 
above the gap. 

The authors wish to acknowledge the ULTRAFAST 
network and the Ministry for Science and Technology of 
Slovenia for supporting part of this work. 



V. APPENDIX 

Analytic solution of Eq. (3) has the form: 

27VLoCexp(-t/r) 



nLit) - 



l-Cexp(-t/T) 



(17) 



Here l/r — 2NlqR + 7. Constant C can be found from 
the following equation (see also Eq.(6)): 



"-l(O) - nL{tr) = rjT^QP- 



(18) 



Combining this two equations one obtains the following 
form for 71^(0): 



nL(0) = iVLo(l+7/2i?iVLo) 



1 



Vngp 



Nl^{l+^/2RNL0?Rtr 



(19) 



This solution reduces to Eq.(8) if 



ynQp 



NlJl+^/2RNLQYRt, 



> 



1 and to Eq.(lO) in the opposite Hmit. 



Figure 1. A photoinduccd transmission signal 
A T I T as a function of time t after photoexcitation in 
YBaaCuaOr-s [T^ = 90K) taken at T = 80 K (points), 
together with the fit (sohd hne). (1) is the basehne sig- 
nal with no pump apphed. (2) is the long-lived signal 
pile-up remaining from previous pulses, (3) (dashed line) 
is the signal due to QP recombination and (4) (dotted 
line) is the long-lived signal remaining after all the QP 
signal has decayed. Signal pile-up contribution and the 
single pulse contribution to the slow photoinduccd signal 
are given by A and ao respectively. 

Figure 2. a) The pump-probe optical diagram. 1 rep- 
resents the pump pulse exciting charge carriers into a 
higher-lying band, 2 the carriers rapidly relax their en- 
ergy to states near the Fermi energy. 3 represents the 
probe pulse, b) A schematic diagram of the terms con- 
tributing to the relaxation of intra-gap states in a super- 
conductor with a gap 2 A. In the case of a CDW gap, the 
e — h pairs take the place of Cooper pairs, so the rate 
equation remains the same. 

Figure 3. a) The temperature dependence of 
the photoinduccd absorption from localized states in 
YBa2Cu306.9 taken from Ref. M (open circles) and Ref. 
[O (solid circles). The solid line is a plot of expres- 
sion (8) with Ac (0) /ksTc = 5 whereas the dashed line 
represents a general solution (Eq. 19) with non-zero 7 
term and additional fourth order term in the expansion 
R in powers of A^ (T) . b) The temperature dependence 
of the photoinduccd reflection in K0.3M0O3 from Ref. p] 
(open circles) compared with the model fit using Eq.(8) 
with Ac{0)/kBTc = 4.8 (solid line). The general solu- 
tion (Eq.(19)) is represented by the dashed line, c) The 
calculated temperature dependence of the photoinduccd 
absorption from localized states in case of T-independent 
pseudogap Eq.(13) with AP/ksTc = 8 and different a/T 
ratios. 
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